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Abstract 

The phenomenon of particle creation within a resonantly vibrating lossy cavity is in- 
vestigated for the example of a massless scalar field at finite temperature. Leakage is 
provided by insertion of a dispersive mirror into a larger ideal cavity. Via the rotat- 
ing wave approximation we demonstrate that for the case of parametric resonance 
the exponential growth of the number of created particles and the strong enhance- 
ment at finite temperatures are preserved in the presence of reasonable losses. The 
relevance for experimental tests of quantum radiation via the dynamical Casimir 
effect is addressed. 

Key words: motion-induced particle creation, canonical approach, finite 
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One of the most impressive manifestations of the nontrivial vacuum structure 
in quantum field theory is the Casimir effect. More than fifty years ago Casimir 
[1] predicted that two conducting parallel plates placed in the vacuum will ex- 
perience an attractive force. This effect has been verified experimentally, see 
e.g. [2]. However, its dynamic counterpart with one or both mirrors moving 
and thus inducing phenomena like the effect of particle creation out of the 
vacuum has not yet been observed experimentally. This striking effect has 
been discussed by many authors (for a review see e.g. [3-6] and references 
therein). It has been shown that under resonance conditions (i.e., when one of 
the boundaries undergoes harmonic oscillations at twice the frequency of one 
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of the eigenmodes of the cavity) the phenomenon of parametric resonance will 
occur. In the case of an ideal cavity (i.e., one with perfectly reflecting mirrors) 
this leads to an exponential growth of the particle occupation numbers of the 
resonance modes [3-6] . In this case a resonantly oscillating boundary is known 
to lead to squeezing of the vacuum state causing the creation of particles in- 
side. 

In view of this result an experimental verification of the dynamical Casimir 
effect in principle appears to be rather simple - provided the cavity is vi- 
brating at the resonance frequency for a sufficiently long time. Of course, this 
point of view is too naive since ideal cavities do not exist. Therefore it is 
essential to include effects of losses. Corresponding investigations have been 
performed for example in [7] based on the conformal invariance of the scalar 
field in 1-1-1 space-time dimensions, see [8]. However, these considerations are 
a priori restricted to 1-1-1 dimensions and cannot simply be generalized to 
higher dimensions. In 3-1-1 dimensions the character of quantum radiation (in 
particular the resonance conditions, see e.g. [3]) differs drastically from the 
1+1 dimensional situation since the spectrum of the eigenfrequencies is not 
equidistant anymore. More realistic cavities with losses were considered in [9]. 
There effects of losses were taken into account by virtue of a (static) master 
equation ansatz and were not derived starting from first principles. 
In addition it is necessary to examine the effects of a detuned external vibra- 
tion frequency [6,9], since in an experiment this frequency will always deviate 
from a desired value [10]. 

Furthermore most papers did not include temperature effects. The canonical 
approach has proven to be quite successful and straightforward and it is also 
capable of including temperature effects. It has been demonstrated in [5] that 
these corrections even enhance the effect of particle production in the case of 
an ideal cavity. However, the canonical approach still lacks a generalization 
for leaky cavities. 

The aim of the current article [10] is to examine the generic properties of the 
dynamical Casimir effect in a non-ideal, resonantly vibrating cavity. For that 
purpose we consider a scalar field inside a dynamical leaky cavity, which serves 
as a suitable model system. One simple way of constructing a leaky system 
is to insert a dispersive mirror into an ideal 3 dimensional cavity thus form- 
ing two leaky cavities coupled to each other. The left cavity is then bounded 
by a perfect mirror ai x — a(t) and a dispersive mirror at x = b. It will be 
considered the leaky cavity, whereas the right one (in addition bounded by a 
perfect mirror at x = c) is understood as the (larger) reservoir, see also Fig. 1. 
Similar - but static - Fabry- Perot systems have been considered frequently in 
quantum optics, see e.g. [11]. 

Note that we are assuming a finite reservoir with a discrete spectrum instead 
of an infinite one leading to a continuum of modes. Since in any experimental 
setup the vibrating cavity will be surrounded by walls, etc., this assumption 
should be justified - cf. the remarks after Eq. (13) - or at least the experiment 
can be designed in this way. 
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Fig. 1. Model of a leaky cavity. A large ideal cavity is split up by a dispersive mirror 
into a lossy cavity and a reservoir. The left (ideal) wall of the cavity is vibrating. 

We consider a massless neutral scalar field coupled to the external potential 
[h — c — kB — 1 throughout] 



£=^ [d,^r; t)] [a^$(r; t)] - V{r; t)^\r; t) , (1) 

where the potential V{r; t) simulates the ideal - inducing Dirichlet boundary 
conditions - and the dispersive mirrors. To keep the calculations simple wc use 
an idealized S-type potential as considered in [12] for the internal dispersive 
mirror 



V{x;t) 



jS{x — b) if a{t) < X < c 
oo otherwise 



Note that when considering photons a dispersive mirror can be realized by a 
thin dielectric slab with a very large dielectric constant. This slab could be 
approximated by a space-dependent permittivity e{x) — l-\-a5{x — h) leading 
to a similar eigenmode equation, see e.g. [11]. However, the generalization to 
real photon creation in dynamic Fabry-Perot cavities is not entirely straight- 
forward (cf. [5]) and remains subject to further elaborations. 
According to [12] the parameter 7 enters the transmission and reflection am- 
plitudes at frequency uj via T — uj/{(jJ + ij) and TZ — —i'y/{u! + i'j). We 
expand the field operator into a complete and orthonormal set of eigenfunc- 
tions ^{r, t) = Qfi{t)ffj_{r; t) satisfying the differential equation {2V{r\ t) — 
A}f^{r;t) = Q?^(t)f^{r]t). A separation of variables leads to an equation ac- 
counting for the x-dependence of the eigenf unctions, which can be solved with 
the ansatz (see also[ll]) 



L^sin[0^(a; — a)] : a{t) < x <b (cavity) 
i?^ sin[Q^(c — x)] : b < X < c (reservoir) 

(3) 



Together with the matching conditions at x = b this leads to a transcendental 
equation for the eigenfrequencies 
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- a) 



cot 



Kic-b) . (4) 



For now we will assume the internal mirror to have a large reflectivity, with 
the parameters 7/^ = being small. Expanding in powers of t]^ one finds 

two sets of frequencies (with a relative shift of O {rj^) compared to those of 
ideal cavities), whose insertion in (3) leads to two separate sets of modes - 
mainly concentrated in the system and the reservoir, respectively. We will 
denote the lowest left-dominated frequency obtained from (4) by fl% and to 
avoid confusion we will also introduce the fundamental perturbation parameter 
rj = rji = ^^2/7- following [10] we shall refer to the different modes as 

left-dominated modes [large in the leaky cavity and small in the reservoir: 
i?^'^^^ — O {rj)] and right-dominated modes [large in the reservoir and small 
in the leaky cavity: L'^^^s^^^ — O (rj)]. Note that when considering mirrors with 
a very high transparency one can always increase the accuracy by including 
higher orders in 7] or even solving the associated equation (4) numerically. 
The total Hamiltonian for the scalar field is obtained from (1) using a Legendre 
transformation. Via insertion of the field operator eigenmode expansion $ it 
can be written as a sum (see e.g. [4]) of a free Hamiltonian 
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and a perturbation Hamiltonian Hj = Hj + Hj with 



H^ = W^Kmi^ Hj^Y^p.^Am^- (6) 

Adopting the nomenclature of [4] we refer to the two parts of the perturbation 
Hamiltonian as the squeezing term Hf and the velocity contribution Hj ■ 
The quantity Ar2^(t) = fi^(t) — (f^^)^ denotes the deviation of the time- 
dependent eigenfrequencies f2^(t) from the unperturbed ones (fi°)^. The inter- 
mode coupling matrix M^y{t) = f cPrf^{r] t)f,^{r; t) is anti-symmetric due to 
the ortho-normality of the modes and the Dirichlet boundary conditions at 
the ideal mirrors, see e.g. [4]. Note that both factors vanish for a stationary 
cavity. For the case of parametric resonance this complicated Hamiltonian 
can be simplified considerably by applying the rotating wave approximation 
(RWA), see e.g. [5,13]. Let us assume the resonance case where during the time 
interval [0, T] the left boundary of the cavity undergoes harmonic oscillations 



a{t) = Oo -l- e{b — Oq) sm{u!t) , (7) 

with a small amplitude < e -C 1 and the external vibration frequency uu. 
With the duration of the vibration being sufficiently long loT » 1 - i.e., after 
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many oscillations - one can approximate the time evolution operator in the 
interaction picture 



U{T, 0) = rexp 



— z 



[Hf{t) + HY{t)) dt 



(8) 



in the following way: 

Each term in the series expansion of (8) can be rewritten as a multiple product 
of Hamiltonians without time-ordering and terms involving commutators with 
Heaviside step functions. The latter terms yield strongly oscillating integrands 
[5,10] and their contribution to (8) is therefore comparably small. Performing 
the time-averaging they will be neglected within the RWA as if one would 
naively neglect the time ordering completely. 

The remaining multi-integrals factorize and can be expanded into powers of 
e < 1 and (uT) > 1. Within the RWA all terms of O (e^(tuT)^) with J > K 
are neglected. (Note that terms with J < K do not occur.) The contributions 
with J — K - i.e., exactly the terms in which the oscillations of the creation 
and annihilation operators ~ exp(inop^) [interaction picture] are compensated 
by the external time dependence [Ai7^(t) ~ sin(cjt) and M^y{t) ~ cos(c<jt)] - 
will be kept. In the resonance case (a; = Vto^) the relevant time integrations 
can be approximated by 



sinM)L^a^Al (9) 



Accordingly, the remaining terms in (8) can be re-summated to yield an ef- 
fective time-evolution operator Uefi{T,Q) — exp{—iH^ffT) with the effective 

interaction Hamiltonian H^gT Jq Hf{t)dt + Jq HJ {t)dt. 
This effective interaction Hamiltonian can be calculated as follows: In terms 
of annihilation and creation operators the squeezing interaction Hamiltonian 
Hf decomposes into a sum of terms like Ar2^(t)aj/^(t)aj/^(t), where AVt^^it) ~ 
2e(Q^°)^ sin(a;t) -|- O (e^) in the resonance case. Due to the trivial time de- 
pendence of the creation and annihilation operators in the interaction picture 
a^,(t) = exp(— ifi^t) the time integration of Hf involves many oscillating 
terms, whose time average is rather small compared to that of constant con- 
tributions. Only the terms where the oscillation of Ar2^(t) compensates the 
oscillation of the operators, i.e., where the squeezing resonance condition (see 
also [5,6,9]) 



u = 2QJ (10) 

holds, will be kept within the RWA. The spectrum of the cavity is assumed 
to be well-separated, i.e., the relative distance of the different energy levels of 
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interest is much larger than O (e). We choose the frequency u to be exactly 
twice the unperturbed frequency of the lowest left-dominated mode (fun- 
damental resonance). 

In the series expansion of (8) this leads with the aid of (9) to a time-averaged 
effective squeezing Hamiltonian 



with the squeezing parameter ^ = efi° 

The same procedure can be applied for the velocity interaction Hamilto- 
nian: Here the coupling matrix M^y{t) factorizes into a time-independent 
geometrical factor = J d?r\df^j^/da)fy and the velocity of the bound- 
ary Mi^,y{t) = m^j^i,a{t) + O (e^). However, here the occurrence of inter-mode 
couplings results in an additional different resonance condition (see also [5,6]) 



(12) 

where we are mainly interested in // = L [10]. In general the above resonance 
condition cannot be fulfilled by a left-dominated mode v - just as for a perfect 
cavity [14]. (The situation is completely different in 1+1 dimensions with Hj 
always contributing due to the equidistant spectrum, see e.g. [3,5].) However, 
in contrast to an ideal cavity it can still be fulfilled by some right-dominated 
mode u — R. Accordingly, one similarly finds with (9) a time-averaged effective 
velocity Hamiltonian 



to ■ 



Hes = iX (aUii - o^ay , (13) 

with X = e^L {^/^y^L + \/^l/^r) mLnib - ao)/4 being the velocity pa- 
rameter of the system. It follows from the characteristics of left- and right- 
dominated modes that the geometry factor becomes small rriL^ = O {rf) in the 
limit of a nearly perfectly reflecting mirror which implies that xl ^ = ^ {v) ^ 
1. The resonance condition (12) can in general be fulfilled by many modes - 
but the main contribution to the particle creation is induced by the squeez- 
ing of the L-mode which couples only to the right-dominated i?-mode in our 
considerations. 

Assuming an arbitrarily large reservoir one would of course obtain many R- 
modes coupling to the fundamental L-mode of interest. However, for a finite 
length of the right cavity c — h and a sufficiently long vibration time T ^ c — h 
the resonance condition will single out a finite number of i?-modes only. Since 
the effect of losses is additive [10] to lowest order in 77 - i.e., for high-quality 
cavities - it is sufficient to consider one i?-mode only. 
Accordingly, the total effective interaction Hamiltonian reads 
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(14) 



We want to calculate the expectation value of particle number operators that 
are explicitly time-independent in the interaction picture: 



{N,{T)) Tr {UI^{T, 0)N,U,^{T, 0)po} , 



(15) 



where po = exp{—(3HQ)/Z denotes the initial statistical operator of the canon- 
ical ensemble. 

Since the Ues is unitary one can introduce new perturbation time (T) depen- 
dent ladder operators 



(16) 



and solve for their T-dependence. This can be done by defining a 4-dimensional 
column vector (cf. [10,6]) 



Via differentiating (16) one can derive an equation for x 



(17) 



dx 
dT 



— I 



(18) 



Since we restrict ourselves to two coupling modes and a quadratic Hamiltonian, 
A is a 4 X 4 number- valued matrix. More sophisticated couplings will simply 
increase the dimension of A. This differential equation is then solved by the 
formal evolution matrix U.{T) 



x{T) = exp (AT) x{0) = U {T)x{0) , 



(19) 



i.e., initial and final ladder operators are related via a Bogoljubov transfor- 
mation. In the present case the effective interaction Hamiltonian (14) implies 
a very simple form of A 






2^ 


X 


2e 





Ox 


-X 











-X 






(20) 
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with Xi = ± — x'^ being its eigenvalues. The time evolution matrix 
U_{T) = exp(AT) is omitted here for brevity but can certainly be calculated, 
e.g. using some computer algebra system. Considering the time evolution of the 
new annihilation and creation operators x{T) — U_ (T)x(O) one finds that the 
expectation values of particle number operators - in particular (Nl) — {x2Xi) 
- can be calculated via inserting x{T) leading to a bilinear form. The full 
response function turns out to be a combination of matrix elements of U.{T) 
[10] 

{Nl{T)) = {UuU21 + UuU23) 

+ {UllU22 + U^2U2l){Nl) 

+ {Ul3U24 + U,^U23){K), (21) 

and similarly for particles created in the reservoir with (Nr) = (x^x^). The 
temperature dependence is taken into account by assuming a Bose-Einstein 
distribution for the initial particle occupation numbers 



of the resonance modes. Taking the limit of f] ^ (ideal cavity) one recovers 
the results found by other authors [3,5,6]. The explicit expression [10] found 
for (iVi(T)) does reflect a purely exponential growth of the particle numbers 
in the resonance modes as long as x < ^, since the time-dependence of U.{T) 
is governed by the eigenvalues of A via exp(AjT'). With x = O {^^lV) this also 
leads to an upper bound for the mirror transmittance r] above which (corre- 
sponding to a highly transparent mirror) one finds oscillations [10]. 
In order to quantify the order of magnitudes let us specify the relevant pa- 
rameters: A cavity with a typical size of A 1 cm would have a fundamental 
resonance frequency of Ql ~ 150 GHz i.e., the coupling right dominated mode 
has a frequency of i^^j = 3^2^ ~ 450 GHz. We will assume a dimensionless vi- 
bration amplitude of e = 10^^, see also [9]. 

Consequently one would have to sustain the vibrations over an interval of 
several milliseconds in order to create a significant number of particles. But 
even after only one millisecond - i.e. ~ 10^ periods - a classical estimate 
based on a quality factor of Q = 10^ [9] would indicate drastic energy losses. 
However, our calculations based on a complete quantum treatment show that 
the effects of losses do not drastically modify the exponential growth of the 
created particles' number as long as 77 <^ 1. Of course our calculations are 
based on the assumption that the larger cavity (including both reservoir and 
the leaky cavity) is perfectly conducting and that the RWA-conditions such 
as T 3> c — 6 hold. The resulting error, however, is of O {Q~'^) and therefore 
- even classically - certainly negligible for 3> 10^, e.g. Q = 10^. Conse- 
quently, the experimental verification of the dynamical Casimir effect could 
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be facilitated by a configTiration where the vibrating cavity is enclosed by a 
larger one. A classical estimate of the quality factor Q for our cavity yields 
Q = 27r/\Tf = 2tt[1 + (7/0;)^] = 0{r]^'^). Accordingly, the assumption of 
Q — 10^, i.e. T) = C(10~^) is completely sufficient to justify our approxima- 
tions. 

Note that a cavity at finite temperatures is even advantageous provided the 
cavity is still nearly ideal for the characteristic thermal wavelength. 
In summary the experimental verification of the dynamical Casimir effect 
could become feasible with lower cavity quahty factors than in [9] (e.g. Q — 
O (10^), which is experimentally achievable) - provided the shift of the reso- 
nance frequency by —77/(26 — 2a) with rj — O (10~^) is taken into account. 

The authors are indebted to A. Calogeracos, V. V. Dodonov, and D. A. R. Dalvit 
for fruitful discussions. R. S. is supported by the Alexander von Humboldt 
foundation and NSERC. Financial support by BMBF and GSI is gratefully 
acknowledged. 
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